3 The effect of dielectric and magnetic media on 
electric and magnetic fields 



3.1 Polarization 

The terrestrial environment is characterized by dielectric media {e.g., air, water) 
which are, for the most part, electrically neutral, since they are made up of neutral 
atoms and molecules. However, if these atoms and molecules are placed in an 
electric field they tend to polarize. Suppose that when a given neutral molecule 
is placed in an electric field E the centre of charge of its constituent electrons 
(whose total charge is —q) is displaced by a distance —r with respect to the centre 
of charge of its constituent atomic nuclei. The dipole moment of the molecule is 
defined p = qr. If there are TV such molecules per unit volume then the electric 
polarization P [i.e., the dipole moment per unit volume) is given by P = Np. 
More generally, 

P{r) = ^Ni{pi), (3.1) 

i 

where {pi) is the average dipole moment of the ith. type of molecule in the vicinity 
of point r, and iVj is the average number of such molecules per unit volume at r. 

It is easily demonstrated that any divergence of the polarization field P{r) 
gives rise to an effective charge density pb in the medium. In fact, 

Pb = -VP. (3.2) 

This charge density is attributable to bound charges {i.e., charges which arise 
from the polarization of neutral atoms), and is usually distinguished from the 
charge density pf due to free charges, which represents a net surplus or deficit of 
electrons in the medium. Thus, the total charge density p in the medium is 

P = Pf+ Pb- (3.3) 

It must be emphasized that both terms in this equation represent real physical 
charge. Nevertheless, it is useful to make the distinction between bound and 
free charges, especially when it comes to working out the energy associated with 
electric fields in dielectric media. 
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Gauss' law takes the differential form 

V.E = ^= Pl±P^. (3.4) 

This expression can be rearranged to give 

V-£> = P/, (3.5) 

where 

D = toE^P (3.6) 

is termed the electric displacement, and has the same dimensions as P (dipole 
moment per unit volume). The divergence theorem tells us that 

(3.7) 




In other words, the flux of D out of some closed surface S is equal to the total 
free charge enclosed within that surface. Unlike the electric fleld E (which is the 
force acting on unit charge) or the polarization P (the dipole moment per unit 
volume), the electric displacement D has no clear physical meaning. The only 
reason for introducing it is that it enables us to calculate fields in the presence 
of dielectric materials without first having to know the distribution of polarized 
charges. However, this is only possible if we have a constitutive relation connecting 
E and D. It is conventional to assume that the induced polarization P is directly 
proportional to the electric field E, so that 

P = eoXeE, (3.8) 

where Xe is termed the electric susceptibility of the medium. It follows that 

D = tQtE, (3.9) 

where 

e=l + Xe (3.10) 

is termed the dielectric constant or relative permittivity of the medium. (Likewise, 
Co is termed the permittivity of free space.) It follows from Eqs. (3.5) and (3.9) 
that 

W-E = ^. (3.11) 
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Thus, the electric fields produced by free charges in a dielectric medium are anal- 
ogous to those produced by the same charges in a vacuum, except that they are 
reduced by a factor e. This reduction can be understood in terms of a polarization 
of the atoms or molecules of the dielectric medium that produces electric fields 
in opposition to that of given charge. One immediate consequence is that the 
capacitance of a capacitor is increased by a factor e if the empty space between 
the electrodes is filled with a dielectric medium of dielectric constant e (assuming 
that fringing fields can be neglected). 

It must be understood that Eqs. (3.8) (3.11) are just an approximation which 
is generally found to hold under terrestrial conditions (provided that the fields are 
not too large) for isotropic media. For anisotropic media {e.g., crystals) Eq. (3.9) 
generalizes to 

D = eoeE, (3.12) 

where e is a second rank tensor known as the dielectric tensor. For strong elec- 
tric fields D ceases to vary linearly with E. Indeed, for sufficiently strong elec- 
tric fields neutral molecules are disrupted and the whole concept of a dielectric 
medium becomes meaningless. 

3.2 Boundary conditions for E and D 

When the space near a set of charges contains dielectric material of non-uniform 
dielectric constant then the electric field no longer has the same form as in vac- 
uum. Suppose, for example, that the space is occupied by two dielectric media 
whose uniform dielectric constants are ei and 62- What are the matching condi- 
tions on E and D at the boundary between the two media? 

Imagine a Gaussian pill-box enclosing part of the boundary surface between 
the two media. The thickness of the pill-box is allowed to tend towards zero, so 
that the only contribution to the outward flux of D comes from its flat faces. 
These faces are parallel to the bounding surface and lie in each of the two media. 
Their outward normals are dSi (in medium 1) and dS2, where dSi = —dS2- 
Assuming that there is no free charge inside the disk (which is reasonable in the 
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limit where the volume of the disk tends towards zero), then Eq. (3.7) yields 



Di-dSi-\- D2-dS2 = 



(3.13) 



where Di is the electric displacement in medium 1 at the boundary with medium 
2, etc. The above equation can be rewritten 



where n2i is the normal to the boundary surface, directed from medium 1 to 
medium 2. If the fields and charges are non time varying then Maxwell's equa- 
tions yield VA-E = 0, which give the familiar boundary condition (obtained by 
integrating around a small loop which straddles the boundary surface) 



In other word, the normal component of the electric displacement and the tan- 
gential component of the electric field are both continuous across any boundary 
between two dielectric materials. 

3.3 Boundary value problems with dielectrics - I 

Consider a point charge q embedded in a semi-infinite dielectric ei a distance d 
away from a plane interface which separates the first medium from another semi- 
infinite dielectric 62- The interface is assumed to coincide with the plane z = 0. 
We need to find solutions to the equations 



(i^2-i^i)-n2i = 0, 



(3.14) 



(E2- El) An2i = 0. 



(3.15) 



P 



(3.16) 



for ^ > 0, 



e2V-E 







(3.17) 



for 2; < 0, and 



VA^ = 



(3.18) 
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everywhere, subject to the boundary conditions at ^ = that 

eiE,{z = 0+) = e2E,{z = 0-), (3.19a) 

E,(z = 0+) = E,{z = 0-), (3.19b) 

Ey{z = 0+) = Ey{z = 0-). (3.19c) 



In order to solve this problem we will employ a slightly modified form of the 
well known method of images. Since V A .E = everywhere, the electric field can 
be written in terms of a scalar potential. So, E = — V0. Consider the region 
^ > 0. Let us assume that the scalar potential in this region is the same as that 
obtained when the whole of space is filled with the dielectric ei and, in addition to 
the real charge q at position A, there is a second charge q' at the image position 
A' (see diagram). If this is the case then the potential at some point P in the 
region ^ > is given by 

> °) = (i + 1 ) • 

where Ri = p'^ + {d — z)'^ and R2 = -^Z p'^ + {d -\- z^, when written in terms of 
cylindrical polar coordinates {p,ip,z). Note that the potential (3.20) clearly is a 
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solution of Eq. (3.16) in the region ^ > 0. It gives V-E = 0, with the appropriate 
singularity at the position of the point charge q. 

Consider the region z < 0. Let us assume that the scalar potential in this 
region is the same as that obtained when the whole of space is filled with the 
dielectric 62 and a charge q" is located at the point A. If this is the case then the 
potential in this region is given by 

The above potential is clearly a solution of Eq. (3.17) in the region ^ < 0. It gives 
V -E = 0, with no singularities. 

It now remains to choose q' and q" in such a manner that the boundary 
conditions (3.19) are satisfied. The boundary conditions (3.19b) and (3.19c) are 
obviously satisfied if the scalar potential is continuous at the interface between 
the two dielectric media: 

(t){z = 0+) = 0(^ = 0-). (3.22) 

The boundary condition (3.19a) implies a jump in the normal derivative of the 
scalar potential across the interface: 

d<l>(z = 0+) d4>(z = 0-) 

— — - — m — • ^^■^^> 



The first matching condition yields 

q^q' q" 



(3.24) 



whereas the second yields 



q-q' = q". (3.25) 



Here, use has been made of 

d f 1\ d f 1 \ d 



(3.26) 
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Equations (3.24) and (3.25) imply that 



9' = -f^^)*' (3.27a) 

«" = f^r-l«- (3.27b) 

V 62 -V t\ ) 



The polarization charge density is given by = — V ■ P, However, inside 
either dielectric P = t^Xe^i so V-P = eoXeV^E = 0, except at the point charge 
q. Thus, there is zero polarization charge density in either dielectric medium. At 
the interface Xe takes a discontinuous jump, 

AXe = ei - €2. (3.28) 

This implies that there is a polarization charge sheet on the interface between 
the two dielectric media. In fact, 

^Tpoi = -(P2-Pi)-n2i, (3.29) 

where 7121 is a unit normal to the interface pointing from medium 1 to medium 
2 (i.e., along the positive 2;-axis). Since 

Pi = eo{ei - 1)E = -eo{ei - 1)V0 (3.30) 

in either medium, it is easy to demonstrate that 

Q €2-61 d 

27rei(e2+ei)(p2 + d2)3/2- ^^'^^^ 

In the limit 62 ^ ei, the dielectric 62 behaves like a conducting medium [i.e., 
E ^ in the region z < 0), and the polarization surface charge density on the 
interface approaches that obtained in the case when the plane z = coincides 
with a conducting surface. 

The above method can clearly be generalized to deal with problems involving 
many point charges in the presence of many different dielectric media whose 
interfaces form parallel planes. 
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3.4 Boundary value problems with dielectrics - II 

Consider a plane slab of dielectric e lying between ^ = and z = b. Suppose that 
this slab is placed in a uniform ^-directed electric field of strength Eq. What is 
the field strength Ei inside the slab? 

Since there are no free charges and this is a one-dimensional problem, it is 
clear from Eq. (3.5) that the electric displacement D is the same in both the 
dielectric slab and the vacuum region which surrounds it. In the vacuum region 
D = gqEq, whereas D = tQtE\ in the dielectric. It follows that 

El = ^. (3.32) 

In other words, the electric field inside the slab is reduced by polarization charges. 
As before, there is zero polarization charge density inside the dielectric. However, 
there is a uniform polarization charge sheet on both surfaces of the slab. It is 
easily demonstrated that 

^^poi(^ = 6) = -cF^o\{z = 0) = eo £^0 (3.33) 
In the limit e 1, the slab acts like a conductor and £^1 — > 0. 
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Let us now generalize this result. Consider a dielectric medium whose dielec- 
tric constant e varies with z. The medium is assumed to be of finite extent and 
is surrounded by a vacuum, so that e(^) — > 1 as \z\ — > oo. Suppose that this 
dielectric is placed in a uniform ^-directed electric field E^. What is the field 
E{z) inside the dielectric? 

We know that the electric displacement inside the dielectric is given by D{z) = 
e^tiz) E{z). We also know from Eq. (3.5) that, since there are no free charges 
and this is a one-dimensional problem. 



dD{z) d\t{z)E{z)\ 



dz ° dz 
Furthermore, E{z) — > Eq as \z\ — > oo. It follows that 

Eq 



0. 



E{z) = 



iz) 



(3.34) 



(3.35) 



Thus, the electric field is inversely proportional to the dielectric constant inside 
the dielectric medium. The polarization charge density inside the dielectric is 
given by 



dE{z) ^ d 
Pb = eo — J — = toEo — 
dz dz 



1 



iz) 



(3.36) 



3.5 Boundary value problems with dielectrics - III 



Suppose that a dielectric sphere of radius a and dielectric constant e is placed in 
a ^-directed electric field of strength Eq (in the absence of the sphere). What is 
the electric field inside and around the sphere? 



Since this is a static problem we can write E = 
charges, so Eqs. (3.5) and (3.9) imply that 

everywhere. The boundary conditions (3.14) and (3.15) reduce to 



dr 



r=a 



dr 



V0. There are no free 

(3.37) 

(3.38a) 



r=a 



+ 



92 



09 



r=a 



09 



(3.38b) 



r=a+ 



Furthermore, 



0(r, 9, cp) -Eq r cos 9 (3.39) 
as r — > 0. Here, (r, 9^ </?) are spherical polar coordinates centred on the sphere. 

Let us search for an axisymmetric solution, = (f){r,9). Since the solutions 
to Poisson's equation are unique, we know that if we can find such a solution 
which satisfies all of the boundary conditions then we can be sure that this is the 
correct solution. Equation (3.37) reduces to 



— sine/ — = 0. 



sm 



9 09 



09 



(3.40) 



Straightforward separation of the variables yields 

oo 

<t>{r, 9) = Y,{Ai + Bi r-('+i))Pz(cos 9), 



(3.41) 



where / is a non- negative integer, the Ai and Bi are arbitrary constants, and Pi{x) 
is a solution to Legendre's equation. 



d 

dx 



(1 - x') 



dx 



+ /(/ + l)P/ = 0, 



(3.42) 



which is single- valued, finite, and continuous in the interval — l<a:;<-|-l. It 
can be demonstrated that Eq. (3.42) only possesses such solutions when I takes 
an integer value. The Pi{x) are known as Legendre polynomials (since they are 
polynomials of order I in x), and are specified by Rodrigues ' formula 



(3.43) 



Since Eq. (3.42) is a Sturm-Liouville type equation, and the Legendre polynomi- 
als satisfy Sturm-Liouville type boundary conditions at a; = ±1, it immediately 
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follows that the Pi (cos 9) are orthogonal functions which form a complete set in 
^-space. The orthogonality relation can be written 

.1 

Pi> (x)Pi (x) dx = 5u. . (3.44) 



2/ + 1 



The Legendre polynomials form a complete set of angular functions, and it 
is easily demonstrated that the and the r~*^'"^^^ form a complete set of radial 
functions. It follows that Eq. (3.41), with the Ai and Bi unspecified, represents a 
completely general axisymmetric solution of Eq. (3.37) which is well behaved in 
^-space. We now need to find values of the Ai and Bi which are consistent with 
the boundary conditions. 

Let us divide space into the regions r < a and r > a. In the former region 

oo 

Hr, e) = ^Ai Pi (cos 6) , (3.45) 
/=o 

where we have rejected the r~^'''^^^ radial solutions because they diverge unphys- 
ically as r — > 0. In the latter region 



oo 



0(r,^) = ^(5/r' + C/r-('+i) )Pz(cos^). (3.46) 
/=o 

However, it is clear from the boundary condition (3.39), and Eq. (3.43), that the 
only non-vanishing Bi is Bi = —Eq. This follows since Pi (cos ^) = cos^. The 
boundary condition (3.38b) (which integrates to give 0(r = a~) = (j){r = a^) for 
a potential which is well behaved in ^-space) gives 

^i = -^o + ^, (3.47) 

and 

^' = ^ (3-48) 

for 1^1. Note that it is appropriate to match the coefficients of the P/(cos^) 
since these functions are orthogonal. The boundary condition (3.38a) yields 

e^i = (3.49) 
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and 

Ci 



tlAi = -(l + l)-^ (3.50) 



a 



for / 7^ 1. Equations (3.48) and (3.50) give Ai = Ci = Qiovl^ 1. Equations (3.47) 
and (3.49) reduce to 

^1 = - ^0' (3-51a) 

Ci = f a^^^o- (3.51b) 



The solution is therefore 



e + 2 



= - ( ^1 ^orcos^ (3.52) 



for r < a, and 

^ e - 1 \ „ 



= -Eorcos6'+ Eo^cos6' (3.53) 

V e + 2 / 



for r > a. 



Equation (3.52) is the potential of a uniform 2;-directed electric field of strength 

Ei = -^Eo. (3.54) 

Note that Ei < Eq provided that e > 1. Thus, the electric field strength is re- 
duced inside the dielectric sphere due to partial shielding by polarization charges. 
Outside the sphere the potential is equivalent to that of the applied field Eq plus 
the field of a point electric dipole, located at the origin and pointing in the z- 
direction, whose dipole moment is 

P = 47reo «'^o- (3.55) 

This dipole moment can be interpreted as the volume integral of the polarization 
P over the sphere. The polarization is 

P = eo(e - 1) El i; = 3eo (7^) ^0 z. (3.56) 
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Since the polarization is uniform there is zero polarization charge density inside 
the sphere. However, there is a polarization charge sheet on the surface of the 
sphere whose density is given by cTpoi = P-r (see Eq. (3.29) ). It follows that 

^^poi = 3eo (^^qr^^ -^0 (3.57) 

The problem of a dielectric cavity of radius a in a dielectric medium with 
dielectric constant e and with an applied electric field Eq parallel to the ^-axis 
can be treated in much the same manner as that of a dielectric sphere. In fact, it 
is easily demonstrated that the results for the cavity can be obtained from those 
for the sphere by making the transformation e ^ 1/e. Thus, the field inside the 
cavity is uniform, parallel to the ^-axis, and of magnitude 

E, = ^^Eo. (3.58) 

Note that Ei > Eq provided that e > 1. The field outside the cavity is the 
original field plus that of a ^-directed dipole, located at the origin, whose dipole 
moment is 

P = -47reo (^^) «'^o- (3.59) 

Here, the negative sign implies that the dipole points in the opposite direction to 
the external field. 



3.6 The energy density within a dielectric medium 

Consider a system of free charges embedded in a dielectric medium. The increase 
in the total energy when a small amount of free charge Spf is added to the system 
is given by 




(3.60) 



where the integral is taken over all space and (j){r) is the electrostatic potential. 
Here, it is assumed that the original charges and the dielectric are held fixed, so 
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that no mechanical work is performed. It follows from Eq. (3.5) that 

SU = J (f)V-SDd^r, (3.61) 

where 6D is the change in the electric displacement associated with the charge 
increment. Now the above equation can also be written 

SU = J V-(4>5D)d^r- j V4>-5Dd^r, (3.62) 



giving 

5U 



j (t)5D-dS- j V4)-5Dd^r, (3.63) 

where use has been made of Gauss's theorem. If the dielectric medium is of finite 
spatial extent then we can neglect the surface term to give 

5U = - J V(f)-5Dd^r = J E-SDd^r. (3.64) 

This energy increment cannot be integrated unless E is a known function of D. 
Let us adopt the conventional approach and assume that D = cqcE, where the 
dielectric constant e is independent of the electric field. The change in energy 
associated with taking the displacement field from zero to D{r) at all points in 
space is given by 

U= 6U= E-SDd^r, (3.65) 

Jo Jo J 



or 




r = \j eoeE^^r, (3.66) 



which reduces to 

U =\j E Dd^r. (3.67) 
Thus, the electrostatic energy density inside a dielectric is given by 

W=^. (3.68) 

This is a standard result which is often quoted in textbooks. Nevertheless, it is 
important to realize that the above formula is only valid in dielectric media in 
which the electric displacement D varies linearly with the electric field E. 
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3.7 The force density within a dielectric medium 



Equation (3.67) was derived by considering a virtual process in which true charges 
are added to a system of charges and dielectrics which are held fixed, so that no 
mechanical work is done against physical displacements. Let us now consider 
a different virtual process in which the physical coordinates of the charges and 
dielectric are given a virtual displacement 6r at each point in space, but no free 
charges are added to the system. Since we are dealing with a conservative system, 
the energy expression (3.67) can still be employed, despite the fact that it was 
derived in terms of another virtual process. The variation in the total electrostatic 
energy 6U when the system undergoes a virtual displacement Sr is related to the 
electrostatic force density / acting within the dielectric medium via 

= - J f-drd^r. (3.69) 

If the medium is moving with a velocity field u then the rate at which electrostatic 
energy is drained from the E and D fields is given by 

dU 



dt 



= - J f ud^r. (3.70) 



Let us now consider the energy increment due to both a change 6pf in the 
free charge distribution and a change 6e in the dielectric constant, caused by the 
displacements. From Eq. (3.67) 

5U=^ f [D'^5{l/e) + 2D-5D/e\ d^r, (3.71) 
2eo J 



or 



6U=-^ j E'^Std^r + j E-6Dd^r. (3.72) 

Here, the first term represents the energy increment due to the change in dielec- 
tric constant associated with the virtual displacements, whereas the second term 
corresponds to the energy increment caused by displacements of the free charges. 
The second term can be written 

J E-5Dd^r = - J V(j)-5Dd^r = j (j)V-5Dd^r = j (j)5pfd^r, (3.73) 
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where surface terms have been neglected. Thus, Eq. (3.72) imphes that 



^^-/ {^^-"-^E'^]d'r. (3.74) 



dt J y dt 2 dt^ 

In order to arrive at an expression for the force density / we need to express 
the time derivatives dp/dt and dt/dt in terms of the velocity field u. This can be 
achieved by adopting a dielectric equation of state; i.e., a relation which gives the 
dependence of the dielectric constant e on the mass density pm- Let us assume 
that t{pm) is a known function. It follows that 

Dt _ dt Dpm , . 

Dt ~ dpm Dt ' ^'^■^^^ 

where 

g.| + «.V (3.76) 

is the total time derivative [i.e., the time derivative in a frame of reference which is 
locally CO- moving with the dielectric.) The hydrodynamic equation of continuity 
of the dielectric is 

^^"^+V-(p^tx)=0, (3.77) 



which implies that 



It follows that 



dt 
Dp, 



Dt 

dt dt 



dt dp 

The conservation equation for the free charges is written 



-PmV-u. (3.78) 



Pm^-u-u-Wt. (3.79) 



dt 

Thus, we can express Eq. (3.74) in the form 

dU 



^^.[pju) = 0. (3.80) 



dt 



d^r. (3.81) 
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Integrating the first term by parts and neglecting any surface contributions, we 
obtain 

(3.82) 



-J 4>^-(pfu) d^r = J pf V4)-ud^r. 



Likewise, 



I dpm 




■ud^r. 



(3.83) 



Thus, Eq. (3.81) becomes 
dU_ I 




■ ud^r. 



(3.84) 



Comparing with Eq. (3.70), we see that the force density inside the dielectric is 
given by 

(3.85) 




The first term in the above equation is the standard electrostatic force density. 
The second term represents a force which appears whenever an inhomogeneous 
dielectric is placed in an electric field. The last term, known as the electro striction 
term, gives a force acting on a dielectric in an inhomogeneous electric field. Note 
that the magnitude of the electrostriction force depends explicitly on the dielectric 
equation of state of the material, through de/dpm- The electrostriction term gives 
zero net force acting on any finite region of dielectric if we can integrate over a 
large enough portion of the dielectric that its extremities lie in a field free region. 
For this reason the term is frequently omitted, since in the calculation of the total 
forces acting on dielectric bodies it usually does not contribute. Note, however, 
that if the electrostriction term is omitted an incorrect pressure variation within 
the dielectric is obtained, even though the total force is given correctly. 



3.8 The Clausius-Mossotti relation 



Let us now investigate what a dielectric equation of state actually looks like. 
Suppose that a dielectric medium is made up of identical molecules which develop 
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a dipole moment 

p = aeoE (3.86) 

when placed in an electric field E. The constant a is called the molecular polar- 
izability. If N is the number density of such molecules then the polarization of 
the medium is 

P = Np = NaeoE, (3.87) 

or 

P = e,E, (3.88) 

where pm is the mass density, Na is Avogadro's number, and M is the molecular 
weight. But, how does the electric field experienced by an individual molecule 
relate to the average electric field in the medium? This is not a trivial question 
since we expect the electric field to vary strongly (on atomic length-scales) inside 
the dielectric. 

Suppose that the dielectric is polarized with a mean electric field Eq which is 
uniform (on macroscopic length-scales) and directed along the ^-axis. Consider 
one of the molecules which constitute the dielectric. Let us draw a sphere of ra- 
dius a about this particular molecule. This is intended to represent the boundary 
between the microscopic and the macroscopic range of phenomena aflfecting the 
molecule. We shall treat the dielectric outside the sphere as a continuous medium 
and the dielectric inside the sphere as a collection of polarized molecules. Ac- 
cording to Eq. (3.29) there is a polarization surface charge of magnitude 

CTpoi = -PcosO (3.89) 

on the inside of the sphere, where (r, 0, (p) are spherical polar coordinates, and 
P = P z = eo(e — 1)Eqz is the uniform polarization of the dielectric. The 
magnitude of at the molecule due to the surface charge is 



where dS = 27ra sin 9 d9 is a surface element of the sphere. It follows that 

E^ = -- I cos^ 9 smOdO = -—. (3.91) 
2eo Jo 3eo 
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It is easily demonstrated that Eq = E^p = at the molecule. Thus, the field at 
the molecule due to the surface charges on the sphere is 

E = ^. (3.92) 



The field due to the individual molecules within the sphere is obtained by 
summing over the dipole fields of these molecules. The electric field at a distance 
r from a dipole p is 



E = 



47reo 



p 3{p-r)r 

/v»3 5 



(3.93) 



It is assumed that the dipole moment of each molecule within the sphere is the 
same, and also that the molecules are evenly distributed throughout the sphere. 
This being the case, the value of E^ at the molecule due to all of the other 
molecules within in the sphere. 



E, 



47reo 



-E 



mols 



P£ _ 3(px xz-\-py yz^-pz z'^) 



(3.94) 



IS zero, smce 



and 



mols 



mols 



mols 



mols 



xy = ^ yz = ^ zx = 0. 



mols 



mols 



mols 



(3.95) 



(3.96) 



It is easily seen that Eq = E^p = 0. Hence, the electric field at the molecule due 
to the other molecules within the sphere vanishes. 

It is clear that the net electric field seen by an individual molecule is 



E — Eq + 



3eo 



(3.97) 



This is larger than the average electric field Eq in the dielectric. The above 
analysis indicates that this effect is ascribable to the long range (rather than the 
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short range) interactions of the molecule with the other molecules in the medium. 
Making use of Eq. (3.88) and the definition P = eo(e — 1)Eq, we obtain 



This is called the Clausius-Mossotti relation. This formula is found to work pretty 
well for a wide class of dielectric liquids and gases. The Clausius-Mossotti relation 
yields 

^ = (i^l)(i±2) 



3.9 Dielectric liquids in electrostatic fields 



Consider the behaviour of an uncharged dielectric liquid placed in an electrostatic 
field. If p is the pressure in the liquid when in equilibrium with the electrostatic 
force density /, then force balance requires that 



Vp=/. 



(3.100) 



It follows from Eq. (3.85) that 
Vp = 



2 ^ 2 \ dpr. 



de 

dprr. 



(3.101) 



We can integrate this equation to give 

rP2 



'Pi 



dp 




^9 de 




^9 de 


,) 


Pm 


~ 2 1^ 


dpm. 


2 


dpm. 





(3.102) 



where 1 and 2 refer to two general points in the liquid. Here, it is assumed that 
the liquid possesses an equation of state, so that p = p{pm)- If the liquid is 
essentially incompressible [pm — constant) then 



P2 - Pi = 



m 



de 

dpm 



(3.103) 
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Finally, if the liquid obeys the Clausius-Mossotti relation then 



P2 -Pl 



eoE^ (6-l)(e + 2)^' 



(3.104) 



According to Eqs. (3.54) and (3.104), if a sphere of dielectric liquid is placed 
in a uniform electric field Eq then the pressure inside the liquid takes the constant 
value 

P=l'oEo''^. (3.105) 

It is clear that the electrostatic forces acting on the dielectric are all concentrated 
at the edge of the sphere and are directed radially inwards; i.e., the dielectric is 
compressed by the external electric field. This is a somewhat surprising result 
since the electrostatic forces acting on a rigid conducting sphere are concentrated 
at the edge of the sphere but are directed radially outwards. We might expect 
these two cases to give the same result in the limit e — > oo. The reason that 
this does not occur is because a dielectric liquid is slightly compressible and is, 
therefore, subject to an electrostriction force. There is no electrostriction force 
for the case of a completely rigid body. In fact, the force density inside a rigid 
dielectric (for which V-it = 0) is given by Eq. (3.85) with the third term (the 
electrostriction term) missing. It is easily seen that the force exerted by an electric 
field on a rigid dielectric is directed outwards and approaches that exerted on a 
rigid conductor in the limit e — > 0. 

As is well known, when a pair of charged (parallel plane) capacitor plates are 
dipped into a dielectric liquid the liquid is drawn up between the plates to some 
extent. Let us examine this effect. We can, without loss of generality, assume 
that the transition from dielectric to vacuum takes place in a continuous manner. 
Consider the electrostatic pressure difference between a point A lying just above 
the surface of the liquid in between the plates and a point B lying just above the 
surface of the liquid well away from the capacitor where E = 0. The pressure 
difference is given by 

PA-PB = - f f dl. (3.106) 

J A 

Note, however, that the Clausius-Mossotti relation yields dejdpm = at both A 
and B, since e = 1 in a vacuum (see Eq. (3.99) ). Thus, it is clear from Eq. (3.85) 
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that the electrostriction term makes no contribution to the hne integral (3.106). 
It follows that 



PA-PB = 



L 



B 



E^Ve-dl. 



(3.107) 



The only contribution to this integral comes from the vacuum/dielectric interface 
in the vicinity of point A (since e is constant inside the liquid, and E = in the 
vicinity of point B). Suppose that the electric field at point A has normal and 
tangential (to the surface) components En and Et, respectively. Making use of 
the boundary conditions that Et and eEn are constant across a vacuum/dielectric 
interface, we obtain 



PA-PB = ^ 



Et\e-l)+e'E^\e)J^^ 



giving 



PA-PB = 



(e - 1) 



TP 2 

2 , 



(3.108) 



(3.109) 



This electrostatic pressure difference can be equated to the hydrostatic pressure 
difference pm g hto determine the height h that the liquid rises between the plates. 
At first sight, the above analysis appears to suggest that the dielectric liquid 
is drawn upward by a surface force acting on the vacuum/dielectric interface 
in the region between the plates. In fact, this is far from being the case. A 
brief examination of Eq. (3.104) shows that this surface force is actually directed 
downwards. According to Eq. (3.85), the force which causes the liquid to rise 
between the plates is a volume force which develops in the region of non-uniform 
electric field at the base of the capacitor, where the field splays out between 
the plates. Thus, although we can determine the height to which the fluid rises 
between the plates without reference to the electrostriction force, it is, somewhat 
paradoxically, this force which is actually responsible for supporting the liquid 
against gravity. 

Let us consider another paradox concerning the electrostatic forces exerted in 
a dielectric medium. Suppose that we have two charges embedded in a uniform 
dielectric e. The electric field generated by each charge is the same as that in 
vacuum, except that it is reduced by a factor e. Therefore, we expect that the 
force exerted by one charge on another is the same as that in vacuum, except 
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that it is also reduced by a factor e. Let us examine how this reduction in 
force comes about. Consider a simple example. Suppose that we take a parallel 
plate capacitor and insert a block of solid dielectric between the plates. Suppose, 
further, that there is a small vacuum gap between the faces of the block and each 
of the capacitor plates. Let ±cr be the surface charge densities on each of the 
capacitor plates, and let icr^ be the polarization charge densities which develop on 
the outer faces of the intervening dielectric block. The two layers of polarization 
charge produce equal and opposite electric fields on each plate, and their effects 
therefore cancel each other. Thus, from the point of view of electrical interaction 
alone there would appear to be no change in the force exerted by one capacitor 
plate on the other when a dielectric slab is placed between them (assuming that 
cr remains constant during this process). That is, the force per unit area (which 
is attractive) remains 

/. = |-- (3.110) 

However, in experiments in which a capacitor is submerged in a dielectric liquid 
the force per unit area exerted by one plate on another is observed to decrease to 




(3.111) 



This apparent paradox can be explained by taking into account the difference 
in liquid pressure in the field filled space between the plates and the field free 
region outside the capacitor. This pressure difference is balanced by internal 
elastic forces in the case of the solid dielectric discussed earlier, but is transmitted 
to the plates in the case of the liquid. We can compute the pressure difference 
between a point A on the inside surface of one of the capacitor plates and a point 
B on the outside surface of the same plate using Eq. (3.107). If we neglect end 
effects then the electric field is normal to the plates in the region between the 
plates and is zero everywhere else. Thus, the only contribution to the line integral 
(3.107) comes from the plate/dielectric interface in the vicinity of point A. Using 
Eq. (3.109), we find that 



PA-PB = 




(3.112) 
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where E is the normal field strength between the plates in the absence of dielectric. 
The sum of this pressure force and the purely electrical force (3.110) yields a net 
attractive force per unit area 




(3.113) 



acting between the plates. Thus, any decrease in the forces exerted by charges 
on one another when they are immersed or embedded in some dielectric medium 
can only be understood in terms of mechanical forces transmitted between these 
charges by the medium itself. 



3.10 Magnetization 

All matter is built up out of atoms, and each atom consists of electrons in mo- 
tion. The currents associated with this motion are termed atomic currents. Each 
atomic current is a tiny closed circuit of atomic dimensions, and may therefore 
be appropriately described as a magnetic dipole. If the atomic currents of a given 
atom all flow in the same plane then the atomic dipole moment is directed normal 
to the plane (in the sense given by the right-hand rule) and its magnitude is the 
product of the total circulating current and the area of the current loop. More 
generally, if j{r) is the atomic current density at the point r then the magnetic 
moment of the atom is 

'm=^J rAjd^r, (3.114) 

where the integral is over the volume of the atom. If there are N such atoms or 
molecules per unit volume then the magnetization M {i.e., the magnetic dipole 
moment per unit volume) is given by M = Nm. More generally, 

M(r) = ^iV,{m,), (3.115) 

i 

where (mj) is the average magnetic dipole moment of the ith type of molecule in 
the vicinity of point r, and iVj is the average number of such molecules per unit 
volume at r. 

Consider a general medium which is made up of molecules which are polar- 
izable and possess a net magnetic moment. It is easily demonstrated that any 
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circulation in the magnetization field M{r) gives rise to an effective current den- 
sity jm in the medium. In fact, 

jm = VAM. (3.116) 

This current density is called the magnetization current density, and is usually 
distinguished from the true current density jt , which represents the convection 
of free charges in the medium. In fact, there is a third type of current called a 
polarization current, which is due to the apparent convection of bound charges. 
It is easily demonstrated that the polarization current density jp is given by 

dP 

i. = ^. (3.117) 

Thus, the total current density j in the medium is given by 

dP 

3=3t^y /\M^—. (3.118) 

It must be emphasized that all terms on the right-hand side of this equation 
represent real physical currents, although only the first term is due to the motion 
of real charges (over more than atomic dimensions). 

The Ampere-Maxwell equation takes the form 

dE 

VAS = /ioj + /ioeo-^, (3.119) 

which can also be written 

dD 

VAB = fiojt+fioy/\M + fJ'O-^, (3.120) 

where use has been made of the definition D = toE -\- P. The above expression 
can be rearranged to give 

VAH = jt + ^, (3.121) 

where 

H= M (3.122) 
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is termed the magnetic intensity, and has the same dimensions as M {i.e., mag- 
netic dipole moment per unit volume). In a steady-state situation, Stokes's the- 
orem tell us that 

(f H dl = [ jfdS. (3.123) 
Jc Js 

In other words, the line integral of H around some closed curve is equal to 
the flux of true current through any surface attached to that curve. Unlike the 
magnetic fleld B (which specifles the force ev A B acting on a charge e moving 
with velocity ?;) or the magnetization M (the magnetic dipole moment per unit 
volume), the magnetic intensity H has no clear physical meaning. The only 
reason for introducing it is that it enables us to calculate fields in the presence of 
magnetic materials without first having to know the distribution of magnetization 
currents. However, this is only possible if we possess a constitutive relation 
connecting B and H. 



3.11 Magnetic susceptibility and permeability 

In a large class of materials there exists an approximately linear relationship 
between M and H. If the material is isotropic then 

M = XmH, (3.124) 

where Xm is called the magnetic susceptibility. If Xm is positive the material is 
called paramagnetic, and the magnetic field is strengthened by the presence of 
the material. If Xm is negative then the material is diamagnetic and the magnetic 
field is weakened in the presence of the material. The magnetic susceptibilities 
of paramagnetic and diamagnetic materials are generally extremely small. A few 
sample values are given in Table 1.^° 

A linear relationship between M and H also implies a linear relationship 
between B and H. In fact, we can write 

B = fiH, (3.125) 

■"^^Data obtained from the Handbook of Chemistry and Physics, Chemical Rubber Company 
Press, Baca Raton, FL 
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Material 


Xrn 


Aluminium 

Copper 

Diamond 

Tungsten 
Hydrogen (1 atm) 
Oxygen (1 atm) 
Nitrogen (1 atm) 


2.3 X 10-^ 
-0.98 X 10-^ 
-2.2 X 10-^ 

6.8 X 10-^ 
-0.21 X 10-8 
209.0 X 10-8 
-0.50 X 10-8 



Table 1: Magnetic susceptibilities of some paramagnetic and diamagnetic mate- 
rials at room temperature 

where 

/^ = /^o(l + Xm) (3.126) 

is termed the magnetic permeability of the material in question. (Likewise, //q 
is termed the permeability of free space.) It is clear from Table 1 that the per- 
meabilities of common diamagnetic and paramagnetic materials do not differ 
substantially from that of free space. In fact, to all intents and purposes the 
magnetic properties of such materials can be safely neglected {i.e., /i = jiq). 

3.12 Ferromagnetism 

There is, however, a third class of magnetic materials called ferromagnetic ma- 
terials. Such materials are characterized by a possible permanent magnetization, 
and generally have a profound effect on magnetic fields {i.e., fi/fio ^ 1). Unfor- 
tunately, ferromagnetic materials do not exhibit a linear dependence between M 
and H or B and H, so that we cannot employ Eqs. (3.124) and (3.125) with 
constant values of Xm and ji. It is still expedient to use Eq. (3.125) as the defi- 
nition of fjL, with fi = fi{H), however this practice can lead to difficulties under 
certain circumstances. The permeability of a ferromagnetic material, as defined 
by Eq. (3.125), can vary through the entire range of possible values from zero 
to infinity, and may be either positive or negative. The most sensible approach 
is to consider each problem involving ferromagnetic materials separately, try to 
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H, amp/m 

Figure 2: Magnetization curve and relative permeability of commercial iron (an- 
nealed) 

determine which region of the B-H diagram is important for the particular case 
in hand, and then make approximations appropriate to this region. 

First, let us consider an unmagnetized sample of ferromagnetic material. If 
the magnetic intensity, which is initially zero, is increased monotonically, then 
the B-H relationship traces out a curve such as that shown in Fig. 2. This is 
called a magnetization curve. It is evident that the permeabilities fi derived from 
the curve (according to the rule fi = B/H) are always positive, and show a wide 
range of values. The maximum permeability occurs at the "knee" of the curve. 
In some materials this maximum permeability is as large as 10^ //q- The reason 
for the knee in the curve is that the magnetization M reaches a maximum value 
in the material, so that 

B = fio{H + M) (3.127) 
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Figure 3: Typical hysteresis loop of a ferromagnetic material 



continues to increase at large H only because of the fioH term. The maximum 
value of M is called the saturation magnetization of the material. 

Next, consider a ferromagnetic sample magnetized by the above procedure. 
If the magnetic intensity H is decreased, the B-H relation does not follow back 
down the curve of Fig. 2, but instead moves along a new curve, shown in Fig. 3, 
to the point R. The magnetization, once established, does not disappear with 
the removal of H. In fact, it takes a reversed magnetic intensity to reduce the 
magnetization to zero. If H continues to build up in the reversed direction, then 
M (and hence B) becomes increasingly negative. Finally, when H increases 
again the operating point follows the lower curve of Fig. 3. Thus, the B-H curve 
for increasing H is quite different to that for decreasing H. This phenomenon is 
known as hysteresis. 

The curve of Fig. 3 is called the hysteresis loop of the material in question. The 
value of B at the point R is called the retentivity or remanence. The magnitude 
of H at the point C is called the coercivity. It is evident that /i is negative in the 
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second and fourth quadrants of the diagram and positive in the first and third 
quadrants. The shape of the hysteresis loop depends not only on the nature of the 
ferromagnetic material but also on the maximum value of H to which the material 
is subjected. However, once this maximum value, -ffmax? becomes sufficient to 
produce saturation in the material the hysteresis loop does not change shape with 
any further increase in -ffmax- 

Ferromagnetic materials are used either to channel magnetic flux {e.g., around 
transformer circuits) or as sources of magnetic field (permanent magnets). For 
use as a permanent magnet, the material is first magnetized by placing it in a 
strong magnetic field. However, once the magnet is removed from the external 
field it is subject to a demagnetizing H. Thus, it is vitally important that a 
permanent magnet should possess both a large remanence and a large coercivity. 
As will become clear later on, it is generally a good idea for the ferromagnetic 
materials used to channel magnetic flux around transformer circuits to possess 
small remanences and small coercivities. 



3.13 Boundary conditions for B and H 

What are the matching conditions for B and H at the boundary between two 
media? The governing equations for a steady state situation are 

V-B = 0, (3.128) 

and 

WAH = jt. (3.129) 

Integrating Eq. (3.128) over a Gaussian pill-box enclosing part of the boundary 
surface between the two media gives 

(B2-B1) -7121 = 0, (3.130) 

where 7121 is the unit normal to this surface directed from medium 1 to medium 
2. Integrating Eq. (3.129) around a small loop which straddles the boundary 
surface yields 

(i?2-i?i) An2i = 0, (3.131) 
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assuming that there is no true current sheet flowing in this surface. In general, 
there is a magnetization current sheet flowing in the boundary surface whose 
density is given by 

J^ = n2iA(M2-Mi), (3.132) 

where Mi is the magnetization in medium 1 at the boundary, etc. It is clear that 
the normal component of the magnetic field and the tangential component of the 
magnetic intensity are both continuous across any boundary between magnetic 
materials. 



3.14 Permanent ferromagnets 

Let us consider the magnetic field generated by a distribution of permanent ferro- 
magnets. Let us suppose that the magnets in question are sufiiciently "hard" that 
their magnetization is essentially independent of the applied field for moderate 
field strengths. Such magnets can be treated as if they contain a fixed, specified 
magnetization iV^(r). 

Let us assume that there are no true currents in the problem, so that jt = 0. 
Let us also assume that we are dealing with a steady state situation. Under these 
circumstances Eq. (3.121) reduces to 

VAH = 0. (3.133) 

It follows that we can write 

H = -V^m, (3.134) 
where 0^ is called the magnetic scalar potential Now, we know that 

V-B = fioV-{H + M) = 0. (3.135) 

Equations (3.134) and (3.135) combine to give 

V^^m = -Pm, (3.136) 

where 

Pm = -V-M. (3.137) 
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Thus, the magnetostatic field H is determined by Poisson's equation. We can 
think of pm as an effective magnetic charge density. Of course, this magnetic 
charge has no physical reality. We have only introduced it in order to make 
the problem of the steady magnetic field generated by a set of permanent mag- 
nets look formally the same as that of the steady electric field generated by a 
distribution of charges. 

The unique solution of Poisson's equation, subject to sensible boundary con- 
ditions at infinity, is well known: 

This yields 

If the magnetization field M{r) is well behaved and localized we can integrate 
by parts to obtain 

^"•W = ^/M(r')-V'(^^)dV. (3.140) 



Now 



so our expression for the magnetic potential can be written 

^".W = -iv./^c;V. (3.142) 

Far from the region of non- vanishing magnetization the potential reduces to 

^".W - A^MciV ^ (3.143) 

where m = J M d^r is the total magnetic moment of the distribution. This 
is the scalar potential of a dipole. Thus, an arbitrary localized distribution of 
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magnetization asymptotically produces a dipole magnetic field whose strength is 
determined by the net magnetic moment of the distribution. 

It is often a good approximation to treat the magnetization field A4"(r) as a 
discontinuous quantity. In other words, A4"(r) is specified throughout the "hard" 
ferromagnets in question, and suddenly falls to zero at the boundaries of these 
magnets. Integrating Eq. (3.137) over a Gaussian pill-box which straddles one 
of these boundaries leads to the conclusion that there is an effective magnetic 
surface charge density, 

am = n-M, (3.144) 

on the surface of the ferromagnets, where M is the surface magnetization, and 
n is a unit outward directed normal to the surface. Under these circumstances 
Eq. (3.139) yields 

47r Jy \r — r \ 47r Jg \r — r \ 

where V represents the volume occupied by the magnets and S is the bounding 
surface to V. Here, dS is an outward directed volume element to S. It is clear 
that Eq. (3.145) consists of a volume integral involving the volume magnetic 
charges pm = — V-iVf and a surface integral involving the surface magnetic charges 
am = n-M. If the magnetization is uniform throughout the volume V then the 
first term in the above expression vanishes and only the surface integral makes a 
contribution. 

We can also write B = V A A in order to satisfy V ■ -B = automatically. It 
follows from Eqs. (3.121) and (3.122) that 

V A ii" = V A {B/fio - M) = 0, (3.146) 

which gives 

V^A = -i^o3m, (3.147) 

since jm = V A M. The unique solution to Eq. (3.147), subject to sensible 
boundary conditions at infinity, is very well known: 
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Thus, 



If the magnetization field is discontinuous it is necessary to add a surface 
integral to the above expression. It is straightforward to show that 

47r Jy \r — r'\ 47r Jg \r — r'\ 

It is clear that the above expression consists of a volume integral involving the 
volume magnetization currents jm = V A M and a surface integral involving 
the surface magnetization currents Jm = M A n (see Eq. (3.132) ). If the mag- 
netization field is uniform throughout V then only the surface integral makes a 
contribution. 



3.15 A uniformly magnetized sphere 

Consider a sphere of radius a, with a uniform permanent magnetization M = 
Mq z, surrounded by a vacuum region. The simplest way of solving this problem 
is in terms of the scalar magnetic potential introduced in Eq. (3.134). From 
Eqs. (3.136) and (3.137), it is clear that 0^ satisfies Laplace's equation, 

V^cf>m = 0, (3.151) 

since there is zero volume magnetic charge density in a vacuum or a uniformly 
magnetized magnetic medium. However, according to Eq. (3.144), there is a 
magnetic surface charge density, 

cTm = r-M = Mocos6', (3.152) 

on the surface of the sphere. One of the matching conditions at the surface of 
the sphere is that the tangential component of H must be continuous. It follows 
from Eq. (3.134) that the scalar magnetic potential must be continuous at r = a, 
so that 

(f)m(r = a+) = (f)m(r = a_). (3.153) 
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Integrating Eq. (3.136) over a Gaussian pill-box straddling the surface of the 
sphere yields 

■ n J. 1 r=a+ 



dr 



= — CTm = —Mq cos 9. 



r=a- 



In other words, the magnetic charge sheet on the surface of the sphere gives rise 
to a discontinuity in the radial gradient of the magnetic scalar potential at r = a. 

The most general axisymmetric solution to Eq. (3.151) which satisfies physical 
boundary conditions at r = a and r = oo is 



oo 



0^(r,^) = ^^zr^Pz (cos ^) 



z=o 



for r < a, and 



1=0 



for r > a. The boundary condition (3.153) yields 

Bi = Ai a^^+i 

for all I. The boundary condition (3.154) gives 

{I + 1) Bi 



-lAia}-^ = -Mo6n 



a 



1+2 



for all /, since P/(cos^) = cos^. It follows that 

Ai = Bi = 

for 1^1, and 

Ai = 



Mn 



B. = 



Mna^ 



(3.155) 



(3.156) 



(3.157) 



(3.158) 



(3.159) 



(3.160a) 
(3.160b) 
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Mr,0) = ^cos^ (3.161) 



Mr:e) = -^z:^cose (3.162) 



Thus, 

Moo^ r 

for r < a, and 

Moa^ a 

for r > a. Since there is a uniqueness theorem associated with Poisson's equation, 
we can be sure that this axisymmetric potential is the only solution to the problem 
which satisfies physical boundary conditions at r = and infinity. 

In the vacuum region outside the sphere 

B = fioH = -iioV(f)m- (3.163) 
It is easily demonstrated from Eq. (3.162) that 



B(r > a) = ^ 



m 3{m-r)r 

^3 ^5 



(3.164) 



where 

m=-7ra^M. (3.165) 

This, of course, is the magnetic field of a magnetic dipole m. Not surprisingly, 
the net dipole moment of the sphere is equal to the integral of the magnetization 
M (which is the dipole moment per unit volume) over the volume of the sphere. 

Inside the sphere we have H = —V(j)m and B = fio{H + M), giving 

M 

H=- — , (3.166) 

and 

B=^fioM. (3.167) 

Thus, both the H and B fields are uniform inside the sphere. Note that the 
magnetic intensity is oppositely directed to the magnetization. In other words, 
the H field acts to demagnetize the sphere. How successful it is at achieving 
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Figure 4: Schematic demagnetization curve for a permanent magnet 

this depends on the shape of the hysteresis curve in the negative H and posi- 
tive B quadrant. This curve is sometimes called the demagnetization curve of 
the magnetic material which makes up the sphere. Figure 4 shows a schematic 
demagnetization curve. The curve is characterized by two quantities: the reten- 
tivity Bn [i.e., the residual magnetic field strength at zero magnetic intensity) 
and the coercivity iiqHc {i.e.., the negative magnetic intensity required to demag- 
netize the material: this quantity is conventionally multiplied by fiQ to give it the 
units of magnetic field strength). The operating point (i.e., the values of B and 
HqH inside the sphere) is obtained from the intersection of the demagnetization 
curve and the curve B = fiH. It is clear from Eqs. (3.166) and (3.167) that 



for a uniformly magnetized sphere in the absence of external fields. The magne- 
tization inside the sphere is easily calculated once the operating point has been 
determined. In fact, Mq = B — /iqH. It is clear from Fig. 4 that for a magnetic 
material to be a good permanent magnet it must possess both a large retentivity 
and a large coercivity. A material with a large retentivity but a small coercivity 



(3.168) 
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is unable to retain a significant magnetization in the absence of a strong external 
magnetizing field. 



3.16 A soft iron sphere in a uniform magnetic field 

The opposite extreme to a "hard" ferromagnetic material, which can maintain a 
large remnant magnetization in the absence of external fields, is a "soft" ferro- 
magnetic material, for which the remnant magnetization is relatively small. Let 
us consider a somewhat idealized situation in which the remnant magnetization 
is negligible. In this situation there is no hysteresis, so the B-H relation for the 
material reduces to 



where //(-B) is a single valued function. The most commonly occurring "soft" 
ferromagnetic material is soft iron [i.e., annealed, low impurity iron). 

Consider a sphere of soft iron placed in an initially uniform external field 
Bq = Bq z. The IjLqH and B fields inside the sphere are most easily obtained 
by taking the solutions (3.166) and (3.167) (which are still valid in this case), 
and superimposing on them the uniform field Bq. We are justified in doing this 
because the equations which govern magnetostatic problems are linear. Thus, 
inside the sphere we have 



B = fi{B)H, 



(3.169) 



B 



2 

Bo + - /ioM. 



(3.170b) 



(3.170a) 



Combining Eqs. (3.169) and (3.170) yields 




(3.171) 



with 




(3.172) 
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where, in general, fi = li{B). Clearly, for a highly permeable material [i.e., 
m/a*o ^ I5 which is certainly the case for soft iron) the magnetic field strength 
inside the sphere is approximately three times that of the externally applied field. 
In other words, the magnetic field is amplified inside the sphere. 

The amplification of the magnetic field by a factor three in the high perme- 
ability limit is specific to a sphere. It can be shown that for elongated objects 
{e.g., rods), aligned along the direction of the external field, the amplification 
factor can be considerably larger than this. 

It is important to realize that the magnetization inside a ferromagnetic mate- 
rial cannot increase without limit. The maximum possible value of M is called the 
saturation magnetization, and is usually denoted Ms- Most ferromagnetic mate- 
rials saturate when they are placed in external magnetic fields whose strengths 
are greater than, or of order, one tesla. Suppose that our soft iron sphere first 
attains the saturation magnetization when the unperturbed external magnetic 
field strength is Bg. It follows from (3.170b) and (3.171) (with fi ':$> fig] that 

B = Bq + 2Bs (3.173) 

inside the sphere, for Bq > Bg- In this case, the field amplification factor is 

B B<, 

^ = 1 + 2^- (3.174) 

-DO -DO 

Thus, for Bq ^ Bs the amplification factor approaches unity. We conclude that 
if a ferromagnetic material is placed in an external field which greatly exceeds 
that required to cause saturation then the material effectively loses its magnetic 
properties, so that // ~ //q- Clearly, it is very important to avoid saturating the 
soft magnets used to channel magnetic flux around transformer circuits. This sets 
an upper limit on the magnetic field strengths which can occur in such circuits. 



3.17 Magnetic shielding 

There are many situations, particularly in experimental physics, where it is de- 
sirable to shield a certain region from magnetic fields. This can be achieved by 
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surrounding the region in question by a material of high permeabihty. It is vitally 
important that a material used as a magnetic shield does not develop a permanent 
magnetization in the presence of external fields, otherwise the material itself may 
become a source of magnetic fields. The most effective commercially available 
magnetic shielding material is called Mumetal, and is an alloy of 5% Copper, 2% 
Chromium, 77% Nickel, and 16% Iron. The maximum permeability of Mumetal 
is about 10^ /io- This material also possesses a particularly low retentivity and 
coercivity. Unfortunately, Mumetal is extremely expensive. Let us investigate 
how much of this material is actually required to shield a given region from an 
external magnetic field. 

Consider a spherical shell of magnetic shielding, made up of material of per- 
meability fi, placed in a formerly uniform magnetic field Bq = Bqz. Suppose 
that the inner radius of the shell is a and the outer radius is b. Since there are 
no free currents in the problem, we can write H = — V0m- Furthermore, since 
B = fiH and V--B = 0, it is clear that the magnetic scalar potential satisfies 
Laplace's equation, V'^(j)m = 0, throughout all space. The boundary conditions 
are that the potential must be well behaved at r = and r — > oo, and also that 
the tangential and the normal components of H and B, respectively, must be 
continuous at r = a and r = b. The boundary conditions on H merely imply that 
the scalar potential 0^ must be continuous at r = a and r = b. The boundary 
conditions on B yield 

d(i)m{r = a-) d(j)m(r = a+) /o i^r: > 

fio = fi , (3.175a) 

fio = fi . (3.175b) 



Let us try the following test solution for the magnetic potential: 

(f)m = ^ r COS ^ + COS ^ (3.176) 

for r > b, 

(f>m= (/5r+^)cos^ (3.177) 
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for b > r > a, and 

^m = 5rcose (3.178) 

for r < a. This potential is certainly a solution of Laplace's equation throughout 
space. It yields the uniform magnetic field Bq as r — )• oo, and satisfies physical 
boundary conditions at r = and infinity. Since there is a uniqueness theorem 
associated with Poisson's equation, we can be certain that this potential is the 
correct solution to the problem provided that the arbitrary constants a, /5, etc. 
can be adjusted in such a manner that the boundary conditions at r = a and 
r = b are also satisfied. 



The continuity of 0^ at r = a and r = b requires that 



and 



The boundary conditions (3.175) yield 



Bo a 
Mo b^ 



a-- 



and 

It follows that 
Mo a = 



Be 



2a 



Mo 



27' 



Mo 



-TT =M 



63 



63 



(3.179) 



(3.180) 



(3.181) 



(3.182) 



(2/i + /io)(M-Mo) 



Mo/5 = - 

M07 = - 
Mo<^ = - 



{2fi + /io)(M + 2mo) - 2(a3/63)(^ _ ^^52 
3(2/i + mo)mo 



(6^ - a^)Bo, 



(2/i + /io)(M + 2mo) - 2(a3/63)(^ _ ^^y J 

3(m - Mo)mo " 

{2fi + /io)(M + 2mo) - 2(a3/63)(^ _ 

9mMo 

(2/i + /io)(M + 2mo) - 2(a3/63)(/i - /io)^ . 



Bo, 

a^Bo, 

Bo. 



(3.183a) 
(3.183b) 

(3.183c) 

(3.183d) 
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Consider the limit of a thin, high permeabihty shell for which b = a -\- d, 
d/a <^ 1, and /i//io ^ 1. In this limit, the field inside the shell is given by 

B.|^^B„. (3.184) 

Thus, if ~ for Mumetal, then we can reduce the magnetic field strength 

inside the shell by almost a factor of 1000 using a shell whose thickness is only 
1/100 th of its radius. Clearly, a little Mumetal goes a long way! Note, however, 
that as the external field strength, Bq, is increased, the Mumetal shell eventually 
saturates, and /i//io gradually falls to unity. Thus, extremely strong magnetic 
fields (typically, Bq > 1 tesla) are hardly shielded at all by Mumetal, or similar 
magnetic materials. 



3.18 Magnetic energy 

Consider an electrical conductor. Suppose that a battery with an electromotive 
field E' is feeding energy into this conductor. The energy is either dissipated as 
heat or is used to generate a magnetic field. Ohm's law inside the conductor gives 

jt = c7(E + E'), (3.185) 

where jt is the true current density, a is the conductivity, and E is the inductive 
electric field. Taking the scalar product with jt, we obtain 

E'-jt = ^-^-E-3t. (3.186) 

(T 

The left-hand side of this equation represents the rate at which the battery does 
work on the conductor. The first term on the right-hand side is the rate of Joule 
heating inside the conductor. We tentatively identify the remaining term with 
the rate at which energy is fed into the magnetic field. If all fields are quasi- 
stationary [i.e., slowly varying) then the displacement current can be neglected, 
and the Ampere-Maxwell equation reduces to V A .ff = j*. Substituting this 
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expression into Eq. (3.186) and integrating over all space, we get 

J E'-{V AH)d^r = j d^r - J E-{V AH)d^r. (3.187) 

The last term can be integrated by parts using the relation 

V-(EAH) =H-(V AE)-E-(V AH). (3.188) 
The divergence theorem plus the Faraday-Maxwell equation yield 

J E-{V AH)d^r = - J H~d^r - J {EAH)-dS. (3.189) 

Since E AH falls off at least as fast as 1/r^ in electrostatic and quasi-stationary 
magnetic fields (1/r^ comes from electric monopole fields, and 1/r^ from magnetic 
dipole fields), the surface integral in the above expression can be neglected. Of 
course, this is not the case for radiation fields, for which E and H fall off like 
1/r. Thus, the constraint of "quasi-stationarity" effectively means that the fields 
vary suflficiently slowly that any radiation fields can be neglected. 

The total power expended by the battery can now be written 

I E'-{VAH)d^r = p^^^d^r + Jn—d^r. (3.190) 

The first term on the right-hand side has already been identified as the energy 
loss rate due to Joule heating. The last term is obviously the rate at which energy 
is fed into the magnetic field. The variation 6U in the magnetic field energy can 
therefore be written 

^U = J H SBd^r. (3.191) 

This result is analogous to the result (3.64) for the variation in the energy of an 
electrostatic field. 

In order to make Eq. (3.191) integrable, we must assume a functional relation- 
ship between H and B. For a medium which magnetizes linearly the integration 
can be carried out in much the same manner as Eq. (3.67), to give 

U =\ j H Bd^r. (3.192) 
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Thus, the magnetostatic energy density inside a hnear magnetic material is given 

by 

W = (3.193) 

Unfortunately, most interesting magnetic materials, such as ferromagnets, exhibit 
a nonlinear relationship between H and B. For such materials, Eq. (3.191) can 
only be integrated between definite states, and the result, in general, depends on 
the past history of the sample. For ferromagnets, the integral of Eq. (3.191) has 
a finite, non-zero value when B is integrated around a complete magnetization 
cycle. This cyclic energy loss is given by 



AU= (p H-dBd^r. (3.194) 




In other words, the energy expended per unit volume when a magnetic material 
is carried through a magnetization cycle is equal to the area of its hysteresis loop 
as plotted in a graph of B against H. Thus, it is particularly important to ensure 
that the magnetic materials used to form transformer cores possess hysteresis 
loops with comparatively small areas, otherwise the transformers are likely to be 
extremely inefficient. 
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